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Abstract 

We investigate the effect of a microwave field on a confined two dimensional electron gas which 
contains an insulating region comparable to the Fermi wavelength. The insulating region causes 
the electron wave function to vanish in that region. We describe the insulating region as a static 
vortex. The vortex carry a flux which is determined by vanishing of the charge density of the 
electronic fluid due to the insulating region. The sign of the vorticity for a hole is opposite to 
the vorticity for adding additional electrons. The vorticity gives rise to non-commuting kinetic 
momenta. The two dimensional electron gas is described as fluid with a density which obeys the 
Fermi-Dirac statistics. The presence of the confinement potential gives rise to vanishing kinetic 
momenta in the vicinity of the classical turning points. As a result the Cartesian coordinate do 
not commute and gives rise to a Hall current which in the presence of a modified Fermi-Surface 
caused by the microwave field results in a rectified voltage. Using a Bosonized formulation of the 
two dimensional gas in the presence of insulating regions allows us to compute the rectified current. 
The proposed theory may explain the experimental results recently reported by J. Zhang et al. 

PACS numbers: 



The topology of the ground state wave function plays a crucial role in determining the 
physical properties of a many-particle system. These properties are revealed through the 
quantization rules. It is known that Fermions and Bosons obey different quantization rules, 
while the quantized Hall conductance [1] and the value of the spin-Hall conductivity are a 
result of non-commuting Cartesian coordinates [2]. Similarly the phenomena of quantum 
pumping observed in one-dimensional electronic systems [3-5] is a result of a space-time cycle 
and can be expressed in the language of non-commuting frequency u = id t and coordinate 
x = idk as shown in ref . [6] . 

Recently the phenomena of rectification current I r (V) = [I(V) + I(—V)]/2 has been 
proposed as a DC response to a low-frequency AC square voltage resulted from a strong 2k p 
scattering in a one dimensional Luttinger liquid [12]. 

In a recent experiment [7], a two-dimensional electron gas (2DEG) GaAs with three insu- 
lating antidots has been considered. A microwave field has been applied, and a DC voltage 
has been measured. The experiment has been performed with and without a magnetic field. 
The major result which occurs in the absence of the magnetic field is a change in sign of the 
rectified voltage when the microwave frequency varies from 1.46 GHz to 17.41 GHz. This 
behavior can be understood as being caused by the antidots, which create obstacles for the 
electrons. 

We report in this letter a proposal for rectification. In section A we present a theory 
which show that rectification can be viewed as a result of non-commuting coordinates. In 
section B we present a qualitative model for rectification, namely the presence of vanishing 
wave function is described by a vortex which induces non-commuting kinetic momenta. The 
sign the vorticity is determined by the vanishing of the electronic density. The electronic 
fluid can be seen as a hard core boson which carry flux, the removal of charge caused by 
the insulating region is equivalent to a decrease of flux with respect the flux of the uniform 
fluid. Including in addition a confining potential we obtain regions where the momentum 
vanishes. The combined effect non-commuting kinetic momenta and confinement gives rise 
to non-commuting cartesian coordinates. In section C we use the Bosonization method to 
construct a quantitative theory which gives rise to a set of equations of motion. Constructing 
an iterative solution of this equations reveals the phenomena of rectifications explained in 
sections A and B. 
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A- Rectifications due to non-commuting coordinates 



Due to the existence of the obstacles, the wave function of the electron vanishes in the 
domain of the obstacles. This will give rise to a change in the wave function , \ip >— > 
|$ >= W(K)\ip > where W(K) is the unitary transformation (induced by the obstacle) and 
the coordinate coordinate representation becomes, r = i-J^ — > r = i-^ + W(K)^U(K) 
[1,2,15]. An interesting situation occurs when the wave function |$ > has zero's [1,2,15] 
or points of degeneracy [16] in the momentum space .This gives rise to non -commuting 
coordinates [1,2,16]. As a result we will have a situation where the the commutator [ri,r 2 ] 
of the coordinates is non zero. 

[niK), r 2 {K)\dK x dK 2 = iVL(K)dK l dK 2 (1) 

Using the one particle hamiltonian h = E(K,r) in the presence of an external electric 
field with the commutators [ri(K),Kj] = i5 it j ,[ri(K) , r 2 (K)] = ifl(K) one obtains [2] the 
Heisenberg equations of motions , 



This equations are identical with the one obtained in ref.[15 ] where E(K, r) is the single 
particle energy being in the semi-classical approximation and E 2 (t) the external electric 
field. As a result of the external electric field E 2 (t) changes the velocity changes according 
to eq.2. Using the interaction picture we find, 

ri = ^[ri(t), -er 2 (t)E 2 (t)\ = ^Q(K)E 2 (t) » ^Q(K)V 2 (t)/L (3) 

V 2 (t) is the voltage caused by the external field E 2 (t). The Fermi Dirac occupation 
function p(K,r) = f FD \E(K,r) — eV 2 (t) — E F \ in the presence of the electric field is used 
to sum over all the single particle states . We obtain the current density J\{r) in the % — 1, 

Jl(r) " 6 / J^ h{R)p{R ^ ~lf ^^ (i?)p(i? ' f) ^ (4) 
The result obtained in the last equation follows directly from the non-commuting coordinates 

— * — * — * 

given by tt(K) ^ 0. The current in eq. 1 depends on p(K, f) = Jf.d. [E(K, r) — eV 2 (t) — E F ], 




FIG. 1: (a) Figure la; (b) Figure lb. 

the Fermi- Dirac occupation function in the presence of the external voltage V 2 (t) ~ E 2 (t)L 
. We expand the non equilibrium density p(K,r) to first order in V 2 (t) we obtain the final 
form of the rectified current. Q(K) has dimensions of a frequency and can be replaced with 
the help of the Larmor's theorem, by an effective magnetic field Q(K) = ^^B e ff(K). This 
allows us to replace eq. 4 by the formula. 

= d ^ fl (K)S(E(K,n - E F) oc (V 2 ( t) m„. 



B-A model for non-commuting coordinates 

We consider a two dimensional electron gas (2DEG) in the presence of a parabolic con- 
fining potential V c (r). The 2DEG contains an insulting region of radius D caused by an 
infinite potential Ui(r) ( in the experiment the insulating region this is caused by three 
antidots) see figure la. The effect of the insulating region of radius D causes the electronic 
wave function \ip(r; R) > to vanish for \ f— R\ < D. The spin of the electrons seems not to 
play any significant role, therefore we approximate the 2DEG by a spinless charge system. 
Such a charged electronic system is equivalent to a hard core charged Boson . For Bosonic 

— * 

wave function has zero's which can be described as a vortex centered at f — R . 

We will show that the following properties are essential in order to have non-commuting 
coordinates. 

— * 

1-The vanishing of the wave function for \ f — R\ < D is described by a vortex localized 
at R. 
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2- The many particles will be described in term of a continuous Lagrange formulation 
[14] f(u,t). Here, f(u,t) is the continuous form of r a (t), where "ct " denotes the particular 
particle, a — 1, 2, . . . , N with a density function po(r), which satisfies N/L 2 = f po{u)d 2 u in 
two dimensions (L 2 is the two dimensional area). The coordinate r(u) and the momentum 
P(u) obey canonical commutation rules, [rj(u), Pj(u')] = iSij5 2 (u — u'). 

2 

3- The parabolic confining parabolic potential V c (r) = 7 ^-f f2 provides the confining length 
L F ,see figures la and lb. 

Using the conditions 1-3 we will show that the non-commuting coordinates emerge: 
-The vanishing of the wave function. 

In the literature it was established that the vanishing of the Bosonic wave function gives 
rise to a multivalued phase and vorticity . See in particular the derivation given in ref. 
[9]. The vortex (the insulating region ) gives rise to non-commuting kinetic momenta , 
[n 1 (r), Il 1 (r')] 7^ where , U = K — d9(f;R) and the phase d$(r; R) is caused by the 
localized vortex [8, 9, 10]. This result is obtained in the following way: 

In the presence of a vortex the single particle operator is parametrize as follows: ip(f; R) = 

\lzA e ie{r,R)^ for < D and = e i6(r,R)^ for > D The field ^ ig & 

— * 

regular hard core boson field and 9(f; R) is a multivalued phase. As a result, the Hamiltonian 
ho = ^K 2 + Ui(r) and the field tjj(f; R) are replaced by the transformed Hamiltonian: 

h = ^(K-d6(r;R)) 2 (5) 
2m 

The momentum K is replaced by the kinetic momentum ,H — K — d0(r; R). The deriva- 

tive of the multivalued phase 9(f;R) determines the vector potential A(f;R) = d9(f;R). 

[n!,n 2 ]^o. 

[UAr), n 2 (r')] = iB(r)6(f- f) w i ^^- 6(r- f) (6) 

where B is an effective magnetic field due to the insulation region, which is defined as 
B(r) = V x A(f;R). The sign of the magnetic field B(r) is determined by the 
vorticity. Following the theory presented in ref.8 (see pages 94-99 and 222-227) 
B(r) has positive vorticity since the electronic density vanishes for the region 
\f — R\ > D creating a hole on background density ( see figure 13.1 page 227 in 
ref.8). 
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For the remaining part of this paper we will replace the delta function by a step function 

— * — * 

E(r, R) which takes the value of one for \ f — R\ < D and zero otherwise. 
-The many particle representation. 

In the presence of hard core Bosons (spinless Fermions) the momenta is replaced by 
h~K — > P(u). The static vortex describes the insulating region and modifies the momentum 
operator , 

n(u)=P(u)-d u #(r(u);R) 

Making use of this continuous formulation, we find a similar result as we have for the 
single particles [9]. 

[Mu),U 2 (u')] = l -^^6\u-u') (7) 
Where E(u, R) = 1 for \u — R\ < D and zero otherwise . 
-The confining potential. 

The last ingredient of our theory is provided by the confining potential and the Fermi 
energy. Due to the confining potential the kinetic momentum has to vanishes for particles 
which have the coordinate close to the classical turning point (see figure 16) \r\ ~ Lp , 
Efermi = ^(|fj ~ Lp) ■ This lead to the following constraint problem for the kinetic 
momentum, 

IIi(|u| faL F )\^ >=0 (8) 

and 



n 2 (|w| ^L F )|^>=0 (9) 

The kinetic momentum IIi(|ti| ~ L F ) and n 2 (|-u| ~ Lp) form a second class constraints 
(according to Dirac's definition [11] the commutator of the constraints has to be non-zero 
) [TIi ( ~ Li?),n 2 (|'u| ~ Lp)] ^ if the region |f| ~ Lp overlaps with the vortex region 

— * — * 

\f— R\ < D. For \u — R\ < D the commutator of the kinetic momenta is given by C = 
given by eq.7. 

We define the matrix [Ci i2 (m, u')]^ 1 = [Ili(-u), Il 2 («')]. Using the function E(u,R) (which 
replaces the delta function ) and the eqs.8,9 we obtain 

id^u')}- 1 w C~ 1 E(u,R)E(u\R)5(\r(u) \ — L F )5(\f(u')\ - L F ) (10) 
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The overlapping conditions are given by the conditions: E(u, R) is equal to one for 

— # 

\u — R\ < D and zero otherwise and S(\f(u) — L F )\ describes the condition of the classical 
turning points. Using eq.10 we find according to Dirac's second class constraints [11] the 
following new commutator [,}v, 



[n(u),r 2 (u')] v = [ n {u),r 2 {rf)\- J du>" J du"[n(u), n^')] (C 1>2 (u", u'"))- 1 [U 2 (u"'),r 2 (u')} 
D 2 E(u, R)E(u>, R)S{\r(u)\ - L F )5(\f(u')\ - L F )] 5{7{u) - r{u')) (11) 



~ i 

For \ f — R\ < D we define a field £l(u; R) trough the equation, 

[D 2 E(u, R)E(u>, R)S{\f(u)\ - L F )S(\f(u')\ - L F )\ = Q(u; R) 

This means that £l(u; R) is approximated by VI oc D 2 for \r — R\ < D . 

Eq. 11 shows that the presence of the momentum , II(m) = P(u) — d u $(r(u); R) 
with the constraints given by eqs. 8,9 gives rise to non-commuting coordinates 
[ri(u),r 2 (u')] v ^0. 

Once we have the result that the coordinate do not commute we can use the analysis 
given in eq 4 (and the result for the current I\ derived with the help of equation 4) to 
compute the rectified current, h oc j[j?-V 2 (t) 

This results can be derived in directly using a modified Bosonization method with a non 
-commuting Kack Moody algebra [16,17]. 

C- The continuous formulation for the 2DEG-A Bosonization approach 

Cl-Bosonization for the 2DEG 

We introduce a continuous formulation for the 2DEG many particles system. We replace 
the single particle Hamiltonian htotai{n, t) by a many electron formulation [14]. We introduce 
a continuous representation, namely f(u, t). Here, r(u, t) is the continuous form of r a (t). The 
coordinate and the momentum obey f(u, t = 0) = u and P(u, t = 0) = K The equilibrium 
Fermi-Dirac density is given by p (u) = f 0? If.dX^K 2 + V c (u) - E F ] 

One of the useful description for many electrons in two dimensions is the Bosonization 
method. We will modify this method [16,17] in order to introduce the effect of the vortex 
field and the confining potential V c {u). 
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C2-The Bosonization method in the absence of the insulating region and con- 
fining potential 

In this section we will present the known [16,17] results for a two dimensional interacting 
metal in the absence of the vortex field and confining potential V c {u). Our starting point 
is the Bosonized form of the 2DEG given in ref. [16,17]. 

Hf.s. = \!<Prj cfr> § \^0tds § ^f^ds'T (s, s>; ? - 7') : 5k n (s, ?) Sk n (s>, ?') : 

Where T (s, s'; ~r — ~r') is the Landau function for the two body interaction [17] and the 
notation :: represents the normal order with respect the Fermi Surface. [17]. According 
to ref.[16,17] the F.S. is described by, k f(s,~F) = k° F (s) + 5 k F (s,~F) .The "normal" 
deformation to the F.S. is given by, 5k\\ (s, ~r) = n (s) ■ 5 k p (s, ~r). "s" is the polar angle 
on the F.S. k (s), and n(s) is the normal to the F.S.. The commutation relations for the 
F.S. are, 

[5k\\(s,~?),5k\\(s',l*)] = i(27r) 2 n(s) ■ V5 2 (n(s) ■ ~r -n(s') ■ P^j 

S (t>° (s) - t° F (o) 

C3- The modification of the Bosonization method in the presence of a confin- 
ing potential 

Following ref. [18] ( see the last term of eq.10 in ref. 18 ) we incorporate into the Bosonic 
hamiltonian the effect of the confining and external potentials. We parametrize the Fermi 
surface in terms of the polar angle s = [0 — 2ir] and the coordinate u .The Fermi surface 
momentum K F (s, u) given by the solution , K F {u) = E F — ^^(-u)^. As a result the the 
FERMI SURFACE ( F.S.) excitations is given by , K F (s, ~u) = K° F {s,u) + 5~~k F (s, Ti) 
.The "normal" deformation to the F.S. is given by, 8k\\ (s, ~u) = n(s,u) ■ 5 k F (s, ~u) and 
n(s, u) is the normal to the F.S. as a function of the polar angle s and real space coordinate 
u. 

Following refs. [16,17,18] we obtain the Bosonized hamiltonian for the many 
particles system in the presence of the potentials , V c (u) and time dependent 
external potential U ext (u,t) . 
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H = J d " u j lK f 2 ny )l [ HlK 2m ^ ^) 2 +Vc(u)5k\\(s, u)+(-e)U ext (u, t)5k {l (s, u)] ds 

(12) 

2 

The new part of in the Bosonic hamiltonian is the presence of the confining V c {u) = r -^-u L 
and external potential U ext {u,t). U ext (u,t) is the external microwave radiation field, 
E 2 (t) = -d 2 U ext (u,t) = E c cos(ut + a(t) and E x {t) = -<9if/ ex '(w, t) = 
The commutation relations for the FERMI SU RF ACE in are given by the Kack Moody 
commutation relation [16,17 ] 

6k\\ (s, it) , <JJfe|| (V, = i (2tt) 2 n(s, u)-V5 2 (n(s, u) ■ u - n(s'u') ■ u'^j 5 (k° f (s, u) - K° F (s', u')^j 

(13) 

C4-The Bosonization method in the presence of the insulating region and 
confining potential 

This problem can be investigated using the hamiltonian given in eq.12 supplemented by 
the constraints conditions imposed by the vanishing density, described by a vortex. Using 
the results given in eq.ll one modifies the commutation relations. This modification can 
be viewed as Dirac's bracket due to second class constraints [11 ]. The commutator [, ] is 
replaced by Dirac 's commutators [,]v- The region of vanishing density is described by the 
function E (u) — 1 for \u — R\ < D and zero otherwise . Using [ri(-u), r 2 (-u')]x? ^ we 
find that the Dirac commutator \6k\\ (s, u) ,6k\\ (s', ~u')] v replaces the commutator given 
in equation 13 



[6k\\ (s, ~u) , 5k\\ (s, ~u')] v = [Sk\\ (s, it') , Sk\\ (s, it')] - 
J d 2 z J dV[^||( S ,^),r 1 (^)]([r 1 (^),r 2 (^0b)"M^(^0,^||(^^')] (14) 

The result given by eq. 14 due to non-commuting coordinates [r 1 (u) : r 2 (u')]v 7^ 0. We 
will compare this result with the one for a magnetic field B , \5k\\ (s, it) , Sk\\ (s 1 , u')] B given 
in ref.17. The commutator for the two dimensional densities in the presence of 
a magnetic field B{u) perpendicular to the 2DEG has been obtained in ref. [17]. 

The modified commutator caused by the magnetic field (see eq. 27 in ref. 17) is: 
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6k\\ (s, u) , 5k\\ (s', u'^j = Sk\\ (s, u) , Sk\\ ^s', u'^j 

^B{u)5\u-u')J-MK\s,u)-K\s\u')} 
ti dt(s) 



(15) 



Where B{u) is the magnetic field and = —sin(s)-^ + cos(s)^ is the derivative in 
the tangential direction perpendicular to the vector n(s) ( the normal to the Fermi surface) 
, n( s ) ■ V = cos^s)^ + sin(s)^ 

Using the analogy between the vortex field an the external magnetic field (eq. 
15 ) we can represent equation 14 in terms of the parameters given in eq.ll. 



5k\\ (s, ~u) , 5k\\ ^s', u'^j = 5k\\(s,u) ,6k\\ (^s 1 ,u'^j — 



D 2 *(.)' — ' « <16) 

To conclude the commutator given in equation 16 allows to investigate the 
physics given in the hamiltonian 12. 

We observe that the Dirac commutator, [5k\\ (s, ~u) ,5k\\ (s', ~u')] v ^ is non-zero for 
s 7^ s'\ The Heisenberg equation of motion will be given by the Dirac bracket. Due to 
the fact that different channels s ^ s' do not commute, the application of an external 
electric field in the i — s$ direction will generate a deformation for the F.S. with s ^ sg-, 

Using this formulation we will compute the rectified current. 

-sin(s)^+cos(s)-^ is the derivative according to the tangential direction 



where 



dt(s) 



which is perpendicular to the vector n(s) with, n(s) ■ V = cos(s)-^ + sin(s)^ 

Using this commutation relations given in equation 16 and the hamiltonian given in 
equation 12 we obtain the equation of motion, 



h[ + } = V F (s,u) 



cos(s) 



d5Us,u;t) 



sims) 



d5ku(s,-u; t) 



duo 



' ok\\{s, u; t ) dt + eh c cos{ut + a{t)) 



1 



sin(s) + cos(s)——E(u, R) 



(17) 



2 m Jo 

We have included in the equation of motion a phenomenological relaxation time r for 
the kinetic momentum. This equation shows that the direct effect caused by the electric 
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field is proportional to sin(s) with the maximum contribution at the polar angles, s — | 
and s = 7r + | where Vjr(s, u) is the Fermi velocity. The effect of the vortex is to generates 
a change in the kinetic momentum perpendicular to the external electric field . This part is 
given by the last term. The last term is restricted to \u — R\ < D and represents the vortex 
contribution .This term is maximum for the polar angles s = and s = ir. The maximum 
effect will be obtained in the region close to the classical turning point where the Fermi 
velocity obeys, V F (s, u) = ^K F (s,u) m 0. 

The current density in the i — 1 direction is given by the polar integration of s , [0 — 2n] . 

J^u) = — [ cos(s) \k° f (s, u)5h\ (s, Ht')] ds (18) 

m J (27r) L J 

We introduce the dimensionless parameter 7 = m ^i D 2 which is a function of ^ and D 
the radius of the insulating region. For values of 7 < 1 we can solve iteratively the equation 
of motion and compute the current. 

In the equation for the kinetic momentum we have included a phenomenological relaxation 
time r. This relaxation time will allow to perform times averages .We only keep single 
harmonics and neglect higher harmonics of the microwave field. 

The iterative solution is given as a series in 7 and the microwave amplitude E c i.e. 
5k\\(s, u; t) = 5k^\s, u; t) + r )5k\^\s, u; t).... 

Solving the equation of motion we determines the evolution of the Fermi surface deforma- 
tion in the presence of the microwave field. We substitute the iterative solution Sk\\(s,u;t) 
obtained from eq.17 into the current density formula given by eq.18. 

In order to provide a physical interpretation of our theory we will use physical parameters 
determined by the experiment. In the experiment the electronic density is n ~ 10 15 m~ 2 this 
corresponds to a Fermi energy of E F ~ O.OleV, equivalent to a temperature of T F ~ 120K 
and a Fermi wavelength of \ F — 0.5 x lCT 7 m. For high mobility GaAs, the typical scattering 
time is r ~ lCT n sec, which corresponds to the mean free path / = v F r. The ratio between 
the mean free path and the Fermi wave length obeys the condition, = = -£r^ — 30. 
Therefore, we can neglect multiple scattering effects. When the thermal length is comparable 
with the size of the system L ~ \hermai — (^") 1//2 ^f, one obtains a ballistic system with 
negligible multiple scattering. 

We describe the confined 2DEG of size L as a system with a parabolic confining potential 
V c (r) = r ^-r i which has a "classical turning point" L F determined by the condition F = 

11 




FIG. 2: 

E F . This condition describes the effective physics of a free electron gas of size L < Lp. 
Demanding that Lp is of the order of the thermal wave length Lp ~ Xthermai determines the 
confining frequency Uq = ^(^);£r- For this condition, we obtain a ballistic regime where 
L< L F ~ 1(T 7 - 10~ 6 m, T ~ 1 - 10K, u ~ 10 10 - 10 u Hz and r ~ 10- n sec. In order to 
be able to observe quantum scattering effects caused by the insulating region of radius "D" , 
we require that the wavelength A^ obeys the condition D > Xp ~ 0.5 x 10 _7 m. 

To leading order in 7 < 1 and in second order in the microwave amplitude E c we compute 
the rectified D.C. voltage V^d.c. i n the % — 1 direction. This rectified voltage is defined as 
V\ t D.c. = h/& (c is conductance in the semi classical approximation determined by the 
transport time which is proportional to the scattering time) . The current l\ given by 
h — x I-l ■M^) ^ 2-u w ^ n L ~ The microwave field is expressed in terms of an R.M.S. 
(effective) voltage Vr.m.s. — E c L/\/2 which allows to define a dimensionless voltage in the 
i — 1 direction i>i,d.c — v^'^ ^ = (x) 2 ^^( < ^)' wnere tan((p) = JZ^a with the function G(y?) 
given in figure 2. 

For 2DEG, we use typical parameters used in the experiment [7], i.e. electronic density 
n 10 15 m -2 with a Fermi energy £jr ~ O.Olef , ojq ~ 10 10 — lO 11 ^^, momentum relaxation 
time t 10 _11 sec. and radius of the insulating region £) > ps 0.5 x 10~ 7 m, with 7 ps 0.7. 
We make a single harmonic approximations (neglect terms which oscillate with frequencies 
2a>o,3a;o ...) We have used figure 3 in ref. [7] to extract the voltage changes as a function 
of the microwave field for a zero magnetic field. Figure 3 in ref. [7] shows clearly a change 
of sign when the microwave varies between 1.46 GHz to 34 GHz and vanishes at 17.41 
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GHz. In figure 2, we have plotted our results given by the formula G(ip) as a function of 
the microwave frequency with the rescaled experimental points (see the three points on our 
theory graph). As shown we find a good agreement of our theory with the experimental 
results once we choose Uq = 17 .41 GHz. For frequencies which obeys < ^- < -p^, we 
find good agreement with the experimental results. However, for low frequencies, our theory 
is inadequate and does not fit the experiment. 

In conclusion, we can say that the origin of the rectification is the emergence 
of the non-commuting Cartesian coordinates and the non-commuting density 
excitations are a result of the vortex field accompanied the classical turning 
caused by the confining potential. Using the modified KacK Moody commutations 
rule for the density excitations we find that excitations with different polar angles s become 
coupled. 

Using this theory we have explained the results of the experiment [7] in a region where 
the magnetic field was zero. 
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